Computational Mutagenesis of Antibody Fragments: Disentangling Side Chains from ΔΔG Predictions.

The development of highly potent antibodies and antibody fragments as binding agents holds significant implications in fields such as biosensing and biotherapeutics. Their binding strength is intricately linked to the arrangement and composition of residues at the binding interface. Computational techniques offer a robust means to predict the three-dimensional structure of these complexes and to assess the affinity changes resulting from mutations. Given the interdependence of structure and affinity prediction, our objective here is to disentangle their roles. We aim to evaluate independently six side-chain reconstruction methods and ten binding affinity estimation techniques. This evaluation was pivotal in predicting affinity alterations due to single mutations, a key step in computational affinity maturation protocols. Our analysis focuses on a data set comprising 27 distinct antibody/hen egg white lysozyme complexes, each with crystal structures and experimentally determined binding affinities. Using six different side-chain reconstruction methods, we transformed each structure into its corresponding mutant via in silico single-point mutations. Subsequently, these structures undergo minimization and molecular dynamics simulation. We therefore estimate ΔΔG values based on the original crystal structure, its energy-minimized form, and the ensuing molecular dynamics trajectories. Our research underscores the critical importance of selecting reliable side-chain reconstruction methods and conducting thorough molecular dynamics simulations to accurately predict the impact of mutations. In summary, our study demonstrates that the integration of conformational sampling and scoring is a potent approach to precisely characterizing mutation processes in single-point mutagenesis protocols and crucial for computational antibody design.

Theoretical analysis of divalent cation effects on aptamer recognition of neurotransmitter targets.

Aptamer-based sensing of small molecules such as dopamine and serotonin in the brain, requires characterization of the specific aptamer sequences in solutions mimicking the in vivo environment with physiological ionic concentrations. In particular, divalent cations (Mg2+ and Ca2+) present in brain fluid, have been shown to affect the conformational dynamics of aptamers upon target recognition. Thus, for biosensors that transduce aptamer structure switching as the signal response, it is critical to interrogate the influence of divalent cations on each unique aptamer sequence. Herein, we demonstrate the potential of molecular dynamics (MD) simulations to predict the behaviour of dopamine and serotonin aptamers on sensor surfaces. The simulations enable molecular-level visualization of aptamer conformational changes that, in some cases, are significantly influenced by divalent cations. The correlations of theoretical simulations with experimental findings validate the potential for MD simulations to predict aptamer-specific behaviors on biosensors.

SiteFerret: Beyond Simple Pocket Identification in Proteins.

We present a novel method for the automatic detection of pockets on protein molecular surfaces. The algorithm is based on an ad hoc hierarchical clustering of virtual probe spheres obtained from the geometrical primitives used by the NanoShaper software to build the solvent-excluded molecular surface. The final ranking of putative pockets is based on the Isolation Forest method, an unsupervised learning approach originally developed for anomaly detection. A detailed importance analysis of pocket features provides insight into which geometrical (clustering) and chemical (amino acidic composition) properties characterize a good binding site. The method also provides a segmentation of pockets into smaller subpockets. We prove that subpockets are a convenient representation to pinpoint the binding site with great precision. SiteFerret is outstanding in its versatility, accurately predicting a wide range of binding sites, from those binding small molecules to those binding peptides, including difficult shallow sites.

GSK-3ß Allosteric Inhibition: A Dead End or a New Pharmacological Frontier?

Most kinase inhibitors are designed to bind to highly homologous ATP-binding sites, which leads to promiscuity and possible off-target effects. Allostery is an alternative approach to pursuing selectivity. However, allostery is difficult to exploit due to the wide variety of underlying mechanisms and the potential involvement of long-range conformational effects that are difficult to pinpoint. GSK-3ß is involved in several pathologies. This critical target has an ATP-binding site that is highly homologous with the orthosteric sites of other kinases. Unsurprisingly, there is also great similarity between the ATP-binding sites of GSK-3ß and its isomer, which is not redundant and thus would benefit from selective inhibition. Allostery would also allow for a moderate and tunable inhibition, which is ideal for GSK-3ß, because this target is involved in multiple pathways, some of which must be preserved. However, despite considerable research efforts, only one allosteric GSK-3ß inhibitor has reached the clinic. Moreover, unlike other kinases, there are no X-ray structures of GSK-3ß in complex with allosteric inhibitors in the PDB data bank. This review aims to summarize the state of the art in allosteric GSK-3ß inhibitor investigations, highlighting the aspects that make this target challenging for an allosteric approach.

Replica-exchange optimization of antibody fragments.

In the framework of the rational design of macromolecules capable of binding to a specific target for biosensing applications, we here further develop an evolutionary protocol designed to optimize the binding affinity of protein binders. In particular we focus on the optimization of the binding portion of small antibody fragments known as nanobodies (or VHH) and choose the hen egg white lysozyme (HEWL) as our target. By implementing a replica exchange scheme for this optimization, we show that an initial hit is not needed and similar solutions can be found by either optimizing an already known anti-HEWL VHH or a randomly selected binder (here a VHH selective towards another macromolecule). While we believe that exhaustive searches of the mutation space are most appropriate when only few key residues have to be optimized, in case a lead binder is not available the proposed evolutionary algorithm should be instead the method of choice.

Arbitrary-Shape Dielectric Particles Interacting in the Linearized Poisson-Boltzmann Framework: An Analytical Treatment.

This work considers the interaction of two dielectric particles of arbitrary shape immersed into a solvent containing a dissociated salt and assuming that the linearized Poisson-Boltzmann equation holds. We establish a new general spherical re-expansion result which relies neither on the conventional condition that particle radii are small with respect to the characteristic separating distance between particles nor on any symmetry assumption. This is instrumental in calculating suitable expansion coefficients for the electrostatic potential inside and outside the objects and in constructing small-parameter asymptotic expansions for the potential, the total electrostatic energy, and forces in ascending order of Debye screening. This generalizes a recent result for the case of dielectric spheres to particles of arbitrary shape and builds for the first time a rigorous (exact at the Debye-Hückel level) analytical theory of electrostatic interactions of such particles at arbitrary distances. Numerical tests confirm that the proposed theory may also become especially useful in developing a new class of grid-free, fast, highly scalable solvers.

Modeling of the Electrostatic Interaction and Catalytic Activity of [NiFe] Hydrogenases on a Planar Electrode.

Hydrogenases are a group of enzymes that have caught the interest of researchers in renewable energies, due to their ability to catalyze the redox reaction of hydrogen. The exploitation of hydrogenases in electrochemical devices requires their immobilization on the surface of suitable electrodes, such as graphite. The orientation of the enzyme on the electrode is important to ensure a good flux of electrons to the catalytic center, through an array of iron-sulfur clusters. Here we present a computational approach to determine the possible orientations of a [NiFe] hydrogenase (PDB 1e3d) on a planar electrode, as a function of pH, salinity, and electrode potential. The calculations are based on the solution of the linearized Poisson-Boltzmann equation, using the PyGBe software. The results reveal that electrostatic interactions do not truly immobilize the enzyme on the surface of the electrode, but there is instead a dynamic equilibrium between different orientations. Nonetheless, after averaging over all thermally accessible orientations, we find significant differences related to the solution's salinity and pH, while the effect of the electrode potential is relatively weak. We also combine models for the protein adsoption-desorption equilibria and for the electron transfer between the proteins and the electrode to arrive at a prediction of the electrode's activity as a function of the enzyme concentration.

Chanalyzer: A Computational Geometry Approach for the Analysis of Protein Channel Shape and Dynamics.

Morphological analysis of protein channels is a key step for a thorough understanding of their biological function and mechanism. In this respect, molecular dynamics (MD) is a very powerful tool, enabling the description of relevant biological events at the atomic level, which might elude experimental observations, and pointing to the molecular determinants thereof. In this work, we present a computational geometry-based approach for the characterization of the shape and dynamics of biological ion channels or pores to be used in combination with MD trajectories. This technique relies on the earliest works of Edelsbrunner and on the NanoShaper software, which makes use of the alpha shape theory to build the solvent-excluded surface of a molecular system in an aqueous solution. In this framework, a channel can be simply defined as a cavity with two entrances on the opposite sides of a molecule. Morphological characterization, which includes identification of the main axis, the corresponding local radius, and the detailed description of the global shape of the cavity, is integrated with a physico-chemical description of the surface facing the pore lumen. Remarkably, the possible existence or temporary appearance of fenestrations from the channel interior towards the outer lipid matrix is also accounted for. As a test case, we applied the present approach to the analysis of an engineered protein channel, the mechanosensitive channel of large conductance.

The Mechanistic Understanding of RAD51 Defibrillation: A Critical Step in BRCA2-Mediated DNA Repair by Homologous Recombination.

The cytotoxic action of anticancer drugs can be potentiated by inhibiting DNA repair mechanisms. RAD51 is a crucial protein for genomic stability due to its critical role in the homologous recombination (HR) pathway. BRCA2 assists RAD51 fibrillation and defibrillation in the cytoplasm and nucleus and assists its nuclear transport. BRC4 is a peptide derived from the fourth BRC repeat of BRCA2, and it lacks the nuclear localization sequence. Here, we used BRC4 to (i) reverse RAD51 fibrillation; (ii) avoid the nuclear transport of RAD51; and (iii) inhibit HR and enhance the efficacy of chemotherapeutic treatments. Specifically, using static and dynamic light scattering, transmission electron microscopy, and microscale thermophoresis, we show that BRC4 eroded RAD51 fibrils from their termini through a "domino" mechanism and yielded monomeric RAD51 with a cumulative nanomolar affinity. Using cellular assays (BxPC-3, pancreatic cancer), we show that a myristoylated BRC4 (designed for a more efficient cell entry) abolished the formation of nuclear RAD51 foci. The present study provides a molecular description of RAD51 defibrillation, an essential step in BRCA2-mediated homologous recombination and DNA repair.

Antibody-Antigen Binding Interface Analysis in the Big Data Era.

Antibodies have become the Swiss Army tool for molecular biology and nanotechnology. Their outstanding ability to specifically recognise molecular antigens allows their use in many different applications from medicine to the industry. Moreover, the improvement of conventional structural biology techniques (e.g., X-ray, NMR) as well as the emergence of new ones (e.g., Cryo-EM), have permitted in the last years a notable increase of resolved antibody-antigen structures. This offers a unique opportunity to perform an exhaustive structural analysis of antibody-antigen interfaces by employing the large amount of data available nowadays. To leverage this factor, different geometric as well as chemical descriptors were evaluated to perform a comprehensive characterization.

Smart Approach for the Design of Highly Selective Aptamer-Based Biosensors.

Aptamers are chemically synthesized single-stranded DNA or RNA oligonucleotides widely used nowadays in sensors and nanoscale devices as highly sensitive biorecognition elements. With proper design, aptamers are able to bind to a specific target molecule with high selectivity. To date, the systematic evolution of ligands by exponential enrichment (SELEX) process is employed to isolate aptamers. Nevertheless, this method requires complex and time-consuming procedures. In silico methods comprising machine learning models have been recently proposed to reduce the time and cost of aptamer design. In this work, we present a new in silico approach allowing the generation of highly sensitive and selective RNA aptamers towards a specific target, here represented by ammonium dissolved in water. By using machine learning and bioinformatics tools, a rational design of aptamers is demonstrated. This "smart" SELEX method is experimentally proved by choosing the best five aptamer candidates obtained from the design process and applying them as functional elements in an electrochemical sensor to detect, as the target molecule, ammonium at different concentrations. We observed that the use of five different aptamers leads to a significant difference in the sensor's response. This can be explained by considering the aptamers' conformational change due to their interaction with the target molecule. We studied these conformational changes using a molecular dynamics simulation and suggested a possible explanation of the experimental observations. Finally, electrochemical measurements exposing the same sensors to different molecules were used to confirm the high selectivity of the designed aptamers. The proposed in silico SELEX approach can potentially reduce the cost and the time needed to identify the aptamers and potentially be applied to any target molecule.

A Fast and Interpretable Deep Learning Approach for Accurate Electrostatics-Driven pKa Predictions in Proteins.

Existing computational methods for estimating pKa values in proteins rely on theoretical approximations and lengthy computations. In this work, we use a data set of 6 million theoretically determined pKa shifts to train deep learning models, which are shown to rival the physics-based predictors. These neural networks managed to infer the electrostatic contributions of different chemical groups and learned the importance of solvent exposure and close interactions, including hydrogen bonds. Although trained only using theoretical data, our pKAI+ model displayed the best accuracy in a test set of ∼750 experimental values. Inference times allow speedups of more than 1000× compared to physics-based methods. By combining speed, accuracy, and a reasonable understanding of the underlying physics, our models provide a game-changing solution for fast estimations of macroscopic pKa values from ensembles of microscopic values as well as for many downstream applications such as molecular docking and constant-pH molecular dynamics simulations.

Pandemic drugs at pandemic speed: infrastructure for accelerating COVID-19 drug discovery with hybrid machine learning- and physics-based simulations on high-performance computers.

The race to meet the challenges of the global pandemic has served as a reminder that the existing drug discovery process is expensive, inefficient and slow. There is a major bottleneck screening the vast number of potential small molecules to shortlist lead compounds for antiviral drug development. New opportunities to accelerate drug discovery lie at the interface between machine learning methods, in this case, developed for linear accelerators, and physics-based methods. The two in silico methods, each have their own advantages and limitations which, interestingly, complement each other. Here, we present an innovative infrastructural development that combines both approaches to accelerate drug discovery. The scale of the potential resulting workflow is such that it is dependent on supercomputing to achieve extremely high throughput. We have demonstrated the viability of this workflow for the study of inhibitors for four COVID-19 target proteins and our ability to perform the required large-scale calculations to identify lead antiviral compounds through repurposing on a variety of supercomputers.

Grouping Hypotheses and an Integrated Approach to Testing and Assessment of Nanomaterials Following Oral Ingestion.

The risk assessment of ingested nanomaterials (NMs) is an important issue. Here we present nine integrated approaches to testing and assessment (IATAs) to group ingested NMs following predefined hypotheses. The IATAs are structured as decision trees and tiered testing strategies for each decision node to support a grouping decision. Implications (e.g., regulatory or precautionary) per group are indicated. IATAs integrate information on durability and biopersistence (dissolution kinetics) to specific hazard endpoints, e.g., inflammation and genotoxicity, which are possibly indicative of toxicity. Based on IATAs, groups of similar nanoforms (NFs) of a NM can be formed, such as very slow dissolving, highly biopersistent and systemically toxic NFs. Reference NMs (ZnO, SiO2 and TiO2) along with related NFs are applied as case studies to testing the oral IATAs. Results based on the Tier 1 level suggest a hierarchy of biodurability and biopersistence of TiO2 > SiO2 > ZnO, and are confirmed by in vivo data (Tier 3 level). Interestingly, our analysis suggests that TiO2 and SiO2 NFs are able to induce both local and systemic toxicity along with microbiota dysbiosis and can be grouped according to the tested fate and hazard descriptors. This supports that the decision nodes of the oral IATAs are suitable for classification and assessment of the toxicity of NFs.

Enhanced Molecular Dynamics Method to Efficiently Increase the Discrimination Capability of Computational Protein-Protein Docking.

Protein-protein docking typically consists of the generation of putative binding conformations, which are subsequently ranked by fast heuristic scoring functions. The simplicity of these functions allows for computational efficiency but has severe repercussions on their discrimination capabilities. In this work, we show the effectiveness of suitable descriptors calculated along short scaled molecular dynamics runs in recognizing the nearest-native bound conformation among a set of putative structures generated by the HADDOCK tool for eight protein-protein systems.

Charged dielectric spheres interacting in electrolytic solution: A linearized Poisson-Boltzmann equation model.

We present an analytical theory of electrostatic interactions of two spherical dielectric particles of arbitrary radii and dielectric constants, immersed into a polarizable ionic solvent (assuming that the linearized Poisson-Boltzmann framework holds) and bearing arbitrary charge distributions expanded in multipolar terms. The presented development entails a novel two-center re-expansion analytical theory that expands upon and improves the existing ones, bypassing the conventional expansions in modified Bessel functions. On this basis, we develop a specific matrix formalism that facilitates the construction of asymptotic expansions in ascending order of Debye screening terms of potential coefficients, which are then employed to find exact closed-form expressions for the total electrostatic energy. In particular, this work allows us to explicitly and precisely quantify the k-screened terms of the potential coefficients and mutual interaction energy. Specific cases of monopolar and dipolar distributions are described in particular detail. Comprehensive numerical examples and tests of series convergence and the relative balance of leading and higher-order terms of the mutual interaction energy are presented depending on the inter-particle distance and particles' radii. The results of this work find application in soft matter modeling and, in particular, in computational biophysics and colloid science, where the availability of increasingly larger experimental structures at the atomic-level resolution makes numerical treatment challenging and calls for more efficient expressions and an increased range of validity.

Molecular Recognition by Gold Nanoparticle-Based Receptors as Defined through Surface Morphology and Pockets Fingerprint.

Ligand shell-protected gold nanoparticles can form nanoreceptors that recognize and bind to specific molecules in solution, with numerous potential innovative applications in science and industry. At this stage, the challenge is to rationally design such nanoreceptors to optimize their performance and boost their further development. Toward this aim, we have developed a new computational tool, Nanotron. This allows the analysis of molecular dynamics simulations of ligand shell-protected nanoparticles to define their exact surface morphology and pocket fingerprints of binding cavities in the coating monolayer. Importantly, from dissecting the well-characterized pairing formed by the guest salicylate molecule and specific host nanoreceptors, our work reveals that guest binding at such nanoreceptors occurs via preformed deep pockets in the host. Upon the interaction with the guest, such pockets undergo an induced-fit-like structural optimization for best host-guest fitting. Our findings and methodological advancement will accelerate the rational design of new-generation nanoreceptors.

Tuning Local Hydration Enables a Deeper Understanding of Protein-Ligand Binding: The PP1-Src Kinase Case.

Water plays a key role in biomolecular recognition and binding. Despite the development of several computational and experimental approaches, it is still challenging to comprehensively characterize water-mediated effects on the binding process. Here, we investigate how water affects the binding of Src kinase to one of its inhibitors, PP1. Src kinase is a target for treating several diseases, including cancer. We use biased molecular dynamics simulations, where the hydration of predetermined regions is tuned at will. This computational technique efficiently accelerates the SRC-PP1 binding simulation and allows us to identify several key and yet unexplored aspects of the solvent's role. This study provides a further perspective on the binding phenomenon, which may advance the current drug design approaches for the development of new kinase inhibitors.

The role of histone tails in nucleosome stability: An electrostatic perspective.

We propose a methodology for the study of protein-DNA electrostatic interactions and apply it to clarify the effect of histone tails in nucleosomes. This method can be used to correlate electrostatic interactions to structural and functional features of protein-DNA systems, and can be combined with coarse-grained representations. In particular, we focus on the electrostatic field and resulting forces acting on the DNA. We investigate the electrostatic origins of effects such as different stages in DNA unwrapping, nucleosome destabilization upon histone tail truncation, and the role of specific arginines and lysines undergoing Post-Translational Modifications. We find that the positioning of the histone tails can oppose the attractive pull of the histone core, locally deform the DNA, and tune DNA unwrapping. Small conformational variations in the often overlooked H2A C-terminal tails had significant electrostatic repercussions near the DNA entry and exit sites. The H2A N-terminal tail exerts attractive electrostatic forces towards the histone core in positions where Polymerase II halts its progress. We validate our results with comparisons to previous experimental and computational observations.